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Spin transport experiments in graphene, a single layer of carbon atoms, indicate spin relaxation 
times that are significantly shorter than the theoretical predictions. We investigate experimentally 
whether these short spin relaxation times are due to extrinsic factors, such as spin relaxation caused 
by low impedance contacts, enhanced spin flip processes at the device edges or the presence of 
an aluminium oxide layer on top of graphene in some samples. Lateral spin valve devices using a 
field effect transistor geometry allowed for the investigation of the spin relaxation as a function of 
the charge density, going continuously from metallic hole to electron conduction (charge densities 
of n ~ lO^^cm"'^) via the Dirac charge neutrality point (n ~ 0). The results are quantitatively 
described by a one dimensional spin diffusion model where the spin relaxation via the contacts 
is taken into account. Spin valve experiments for various injector/detector separations and spin 
precession experiments reveal that the longitudinal (Ti) and the transversal (T2) relaxation times 
are similar. The anisotropy of the spin relaxation times ry and rx, when the spins are injected 
parallel or perpendicular to the graphene plane, indicates that the effective spin orbit fields do not 
lie exclusively in the two dimensional graphene plane. Furthermore, the proportionality between 
the spin relaxation time and the momentum relaxation time indicates that the spin relaxation 
mechanism is of the Elliott- Yafet type. For carrier mobilities of 2-5x10^ cm'^/Vs and for graphene 
flakes of 0.1-2 /xm in width, we found spin relaxation times of the order of 50-200 ps, times which 
appear not to be determined by the extrinsic factors mentioned above. 

PACS numbers: 72.25.-b 



I. INTRODUCTION 

Spintronics is the field of electronics that uses the spin 
of the charge carrier for the transport and manipulation 
of information. Preserving a spin state is fundamentally 
limited by the spin orbit interaction which provides the 
mechanism(s) for spin relaxation. Carbon based mate- 
rial systems promise a weak spin orbit interaction (i.e., 
long spin relaxation times) due to the low atomic num- 
ber of carbon. Spin transport investigations have been 
reported for molecular systems and carbon nanotubes, 
see [l| for a review jand more recently for sin gle or multi- 
layer graphene @, i, i, i, i, 0, S i 0, [HI, 111 [H, • 
Combined with the relatively large carrier mobilities 

[HE 

\vf\ , this determines a long spin relaxation length, 
i.e. a long distance which a spin can travel without loos- 
ing its initial orientation. 

In previous experiments we deduced spin relaxation 
times of 150 ps (spin relaxation lengths of 1.5-2 /xm) at 
room temperature @, value confirmed by more recent 
findings (l3 . [l3 | . These relaxation times are shorter than 
expected by at least one order of magnitude, pointing to 
an extrinsic spin relaxation mechanism rather than to an 
intrinsic one. In this study, we present all electrical spin 
valve and spin precession experiments aiming at deter- 
mining the possible causes for such short spin relaxation 
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times. From the experimental point of view, three major 
aspects are identified as discussed below. For theoretical 
investigations of the possible s pin relaxation mechanisms 
we refer to [ll, E, 0, Ull, HI tH Hi, i25j . 

First, in order to combat the conductivity mismatch 
problem (26j and realize efficient electrical spin injection 
and detection in graphene, a thin aluminium oxide layer 
has been introduced in between the ferrornagnetic Co 
injector/detector and graphene [H, [l^, [H, 0, [s^, [3]| . 
Yet, in some experiments the contact resistances are 
comparable in magnitude with the graphene square re- 
sistance over one spin relaxation length and significant 
spin relaxation occurs via the contacts. We quantify 
this effect by using an one dimensional spin diffusion 
model. Second, in all our previous experiments the whole 
graphene layer was covered by the thin aluminum oxide 
layer, a procedure which simplified the fabrication pro- 
cess. Whether this oxide layer, which needs to be only 
underneath the Co electrodes, induces extra spin relax- 
ation was not known. In order to investigate this issue we 
performed experiments where the aluminium oxide barri- 
ers were present only underneath the Co electrodes. The 
third aspect is whether the low spin relaxation times are 
due to spin scattering at the edges of the graphene fiakes. 
Therefore, we performed spin transport experiments us- 
ing narrow structures (down to 100 nm in width, 15-20 
times narrower than the spin relaxation length). 

This paper is organized as follows. In Section [ill 
after a short introduction regarding the functioning of 
a four terminal spin valve device, we present fabrica- 
tion/measurement details and general electrical charac- 
teristics of the devices. In Section lllll we discuss modeling 
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and experiments pertaining the spin valve and spin pre- 
cession experiments. Details of the modeling of spin re- 
laxation via the contacts are presented in the Appendix. 



II. EXPERIMENTS 

A. Graphene lateral spin valve devices 

In Fig. [T] we show scanning electron micrographs 
(SEM) of two graphene flakes contacted by four ferro- 
magnetic Co electrodes (25-40 nm in thickness and differ- 
ent widths of up to 800 nm) . All devices were fabricated 
on an oxidized silicon substrate and have an aluminium 
oxide layer all over the graphene (type I devices) or only 
underneath the Co electrodes (type II and type III de- 
vices). The fabrication procedure which defines the de- 
vice types is described in the next section. 




FIG. 1: SEM of two graphene spin valves (type III devices, 
fabricated by etching) . 

In Fig. [2K) we show a schematic cross section of a lat- 
eral graphene type I spin valve device in the non-local 
measurement geometry (ssl . |33 | . For simplicity, we con- 
sider the outer electrodes (Fl and F4) as being non- 
magnetic. The current set between F2 and Fl creates 
a spin accumulation (imbalance) in graphene at the posi- 
tion of F2, i.e. the electrochemical potentials for the spin 
up and spin down (/i-f,/ij^) channels split [3^ . Fig. [2]d). 
The spin accumulation fis = fJ-i ~ fJ-i diffuses away from 
the injection point {x — 0) and, in the linear regime, 
obeys the Bloch equation 



r 



a; X /x^ = 0, 



(1) 



where D and r represent the spin diffusion constant and 
the spin relaxation time, respectively. The term u: x fi^ 
describes the precession of the spin accumulation in an 
external magnetic field B with the Larmor frequency uj = 
qHbB, with g = 2 the gyromagnetic factor and fiB the 
Bohr magneton. In Eq.[T]we used a single relaxation time 
T for reasons discussed later (see also jsl. Issjp. 

The spin accumulation is probed non-locally by the 
voltage difference Vni between F3 and F4, with F3 placed 
at a distance x ~ L oi the order of A with respect to the 
injection point. Depending on its magnetization orienta- 
tion (parallel/antiparallel) with respect to the spin accu- 
mulation, electrode F3 is sensitive to the electrochemical 
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FIG. 2: a) Schematic cross-section of a type I graphene spin 
valve device measured in the non-local geometry, b) Repre- 
sentation of the electrical spin injection/detection process, c) 
Spin valve measurement (device I. A) showing a 50 SI spin sig- 
nal. Note the sharp switching and the symmetry of the spin 
signal around zero non-local resistance. 



potential of either the spin up channel or the spin down 
one, Fig.[2b). The spin signal is defined as the non-local 
resistance i?„; — Vni/I, where / represents the injection 
current. In the absence of precession a; x /x^ = the 
spin accumulation decays exponentially with the distance 
with the characteristic length A = V Dt, the spin relax- 
ation length. Under the assumption of high impedance 
ferromagnets one can show that 



P^RsqX . L 

Rni — ± — exp 



(2) 



where Eq. 1 from [3J| has been adapted for the two di- 
mensional graphene. Rsq represents the graphene square 
(sheet) resistance, W is the width of the flake and P 
is the spin inject ion/ detect ion efficiencies of the ferro- 
magnetic electrodes. The -l-(— ) signs correspond to the 
parallel (antiparallel) orientation of the magnetization of 
the injector (F2) and the detector (F3). The spin valve 
measurement, see Fig. [2t), consists in monitoring the 
non-local resistance as a function of an external magnetic 
field which manipulates the relative orientation of the in- 
jector/detector magnetizations. The electrodes are engi- 
neered to have different widths so that, due to the shape 
anisotropy, they switch their magnetization orientation 
at different magnetic fields applied along their easy axis 
(the length). The spin valve signal Ai?„; = 2|i?„i| repre- 
sents the change in resistance when the injector/detector 
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magnetization configuration changes from parallel to an- 
tiparallel. The exponential dependence of the spin valve 
signal Ai?„; on the injector-detector separation L allows 
for the extraction of the spin relaxation length, the pa- 
rameter of most interest. 

In the non-local geometry, due to the separation of 
the injection and detection circuits, magnetoresistance 
effects related to the charged current flow, such as 
anisotropic magnetoresistance and Hall effects, do not 
superimpose on the spin signal. Theoretically, the non- 
local resitance is determined only by the spin accumula- 
tion injected at a; = which diffused towards the detector 
(x = L), i.e., there are no other resistances (voltages) on 
top of the non-local resistance (voltage). In most sam- 
ples, however, a small background (frequency and gate 
voltage dependent) is measured. In a local measurement 
geometry, the injection and detection circuits share the 
same pair of electrodes. The spin signals are superim- 
posed on the graphene resistance (typically few kJl) and 
are more difficult to detect, especially if the contact or 
graphene resistance fluctuates. Nevertheless, there are a 
few local measurements reported [1, . 



B. Device fabrication and electrical characteristics 

The devices were fabricated on a dry oxidized, highly 
doped (p=lmf2cm) silicon substrate commercially avail- 
able. The oxide was 500 nm or 300 nm thick. The ap- 
plication of a back gate voltage allows for the control 
of the charge carrier type and density in the graphene 
flakes. The gate electrode was fabricated by etching 
the back side oxide of the wafer followed by the depo- 
sition of a 100 nm thick Ti/Au layer. On this substrate 
we defined a set of Ti/Au markers using electron beam 
lithography (EBL). Next, the substrate was cleaned by 
oxygen plasma in order to remove the polymeric residu- 
als from the EBL step. Then, highly oriented pyrolytic 
graphite (HOPG from GE Advanced Ceramics and SPI) 
was cleaved with a scotch tape and firmly pressed against 
the substrate. The tape is pealed off the substrate leav- 
ing a random distribution of graphite pieces among which 
there are also single layer graphene flakes. Using an op- 
tical microscope, the lowest contrast flakes were selected 
as potential single layer candidates. After that, the sin- 
gle layers were identified using atomic force microscope 
(AFM) measurements. At the same time their position 
with respect to the known position of the Ti/Au markers 
is noted which later on allowed for a precise positioning 
of the ferromagnetic electrodes. In a few cases Raman 
spectroscopy was performed in order to confirm the single 



layer selection performed by AFM [37[. The fabrication 
procedure continued differently for the different types of 
devices as explained below. 

For the type I samples, we continued with the thermal 
evaporation of a thin layer of Al (6 A) in an ultra high 
vacuum (UHV) system with the substrate being liquid 
nitrogen cooled. After that, the Al layer was oxidized for 



minimum 30 minutes in the load lock of the system in 
an oxygen atmosphere of about 100 mbar. Next, ferro- 
magnetic electrodes of different widths were defined using 
EBL followed by electron beam evaporation of Co in a 
high vacuum system. After lift off in hot acetone for typ- 
ically 10 minutes, the sample was glued on a chip carrier 
and electrical connections between the sample and the 
chip carrier were made using ultrasonic wire bonding. 
The sample was then placed in a vacuum container and 
measured. For the type I devices, the graphene flakes 
are completely covered by the thin aluminium oxide. In 
Fig. [3^)-b) we show AFM measurements of a graphene 
flake after the aluminium deposition. In spite of the cryo- 
genic evaporation of Al, the oxide layer on graphene is 
granular in nature. SEM images, not shown here, sup- 
port this observation. 




3 ^jm 



FIG. 3: AFM measurements of two graphene flakes, a) and 
b) height and phase data (600 x 600 nm^) after Al deposition 
(type I device), c) and d) height data (3x3 ^irr?) after etching 
and after annealing in Ar-|-H2 atmosphere (type III device). 

In order to determine the influence of the AlOx layer 
covering graphene on the spin transport properties, de- 
vices where the AlOx layer was deposited only under- 
neath the Co electrodes were fabricated (type II de- 
vices). In this case an EBL step was performed before 
the sequential deposition/oxidation of Al and the ther- 
mal evaporation of Co in the same UHV system. Due to 
the configuration of the UHV system the Al deposition 
on the cooled stage was done under an angle of 30° with 
respect to the sample normal. The deposition was set 
to take place along the electrode length with a precision 
of about 1-2°. This could potentially result in shadow- 
ing effects due to the resist walls, yielding non-uniform 
deposition of Al in the patterned areas. 

The type HI devices were fabricated in order to inves- 
tigate spin transport in narrow graphene ribbons. In this 
case all samples have the AlOx layer only underneath the 
Co contacts (same as type II). However, naturally narrow 
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flakes are difficult to spot under an optical microscope. 
To fabricate devices with widths down to 100 nm, some- 
what larger flakes were selected and then oxygen plasma 
etched, prior to the definition of the composite AlOx/Co 
electrodes. The oxygen plasma etching step often re- 
sulted in a contaminated surface. Even though after 
etching the flakes were annealed in an Ar(95)%:H2(5%) 
atmosphere at 350° for about two hours, large particles 
of unknown chemical composition still remain on the sur- 
face, see Fig. [3J:)-d). 

The electrical measurements were performed using a 
standard AC lock-in technique (/ <17 Hz) in the four 
terminal non-local geometry. The AC current (1-20 fiA) 
was set between a pair of electrodes and the in phase AC 
non-local voltage was measured between the other pair 
of electrodes. The room temperature measurements were 
performed with the sample inserted in a vacuum con- 
tainer, base pressure of ^ 5 x 10~® mbar (measured close 
to the pump), placed in between the poles of a room tem- 
perature electromagnet. The measurements were done in 
vacuum in order to reduce the hysteresis of the graphene 
resistance as a function of the gate voltage and to avoid 
the degradation of the spin signal as a function of time. 
In the case of atmospheric measurement conditions, the 
reduction of the spin signal in time (on a time scale of 
about 24 hours) was accompanied by an increase of the 
contact resistances. The low temperature (4.2 K) mea- 
surements were performed in an Oxford cryostat with the 
sample placed in a dip-stick and using a superconducting 
magnet. 

In Fig. m we compare electrical characteristics of 
graphene and contacts for two devices, one type I and one 
type II. As a function of the gate voltage, the four termi- 
nal measurements of the graphene resistivity (top panels) 
show a maximum which corresponds to the Dirac neu- 
trality point where the average charge density is equal to 
zero. At negative voltages with respect to the Dirac point 
position the conduction takes place via holes, whereas for 
positive voltages the conduction happens via electrons. 
The mobility of all samples presented in this study was in 
the range 2-5x10^ cm^/Vs. A clear correlation between 
the carrier mobilities of the type I (completely covered 
by aluminium oxide) and type II samples (oxide barriers 
only underneath the electrodes) could not be established. 
For the type III devices the carrier mobilities were the 
lowest (2-2.5x10^ cmVVs). 

In the lower panels of Fig. [4] we show three-terminal 
contact resistance measurements as a function of contact 
area, gate voltage and DC current bias. In general, the 
contact resistances of type I and type III devices show a 
large spread in values and no clear scaling with the con- 
tact area. For the type II samples, the spread is less and 
there seems to be some scaling of the contact resistance 
with the area. For all types of contacts, a small and ir- 
regular gate voltage dependence of the contact resistance 
is recorded. Some low resistance contacts show a similar 
behavior to graphene, i.e. the contact resistance is higher 
when the graphene flake is at the charge neutrality point. 
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FIG. 4: Graphene and contact electrical characteristics (dif- 
ferential resistance dV/dl) for type I (left column) and type 
II (right column) devices. 



Most likely, this indicates that part of the graphene layer 
contributes to the contact resistance. With respect to the 
DC bias current dependence, some high resistance con- 
tacts show high values in zero bias conditions whereas 
low resistance contacts show very weak or no DC bias 
dependence. Unfortunately, the large spread in the over- 
all behavior does not allow for a quantitative analysis. 
We take these observations and the AFM measurements 
(Fig. [3]) as evidence that the contacts characteristics are 
determined by a random distribution of regions with in- 
creased transparency (thinner oxide layer or pinholes) in 
contrast to the ideal case of tunnel barriers. 



III. MODELING AND EXPERIMENTS 

A. Modeling the contact resistance effects in the 
spin valve measurements 

Analyzing Fig.l^h-) we see that the spins injected by F2 
do not only diffuse symmetrically to the left and to the 
right in graphene, but can also return into the ferromag- 
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netic electrode where they loose their initial orientation 
very fast (due to the strong spin orbit coupling in the 
ferromagnet) . Therefore, the magnitude of the spin ac- 
cumulation created in graphene is drastically reduced if 
the contact resistance is much lower than the graphene 
resistance over one spin relaxation length. This effect is 
known as the conductivity mismatch problem [2(|]. Tak- 
ing into account the spin current drawn by the electrode, 
see the Appendix for details, one can show that 



± 



2P^RsqX (i?/A)2cxp(-L/A) 



W (l + 2i?/A)2-cxp(-2L/A)' 
where the parameter R given by 

R^^W 



(3) 



(4) 



represents spin relaxation due to the finite contact re- 
sistances Rc of the injector (F2) and detector (F3). In 
fact, R/X represents the ratio between the contact re- 
sistance and the graphene resistance over one spin relax- 
ation length. Eq.[3]is in agreement with [s^l given our as- 
sumptions of negligible ferromagnet resistances with re- 
spect to the contact and graphene resistances and small 
polarizations P of the injector/detector electrodes. 

In the limit of high impedance contacts {Rc — > oo) 
Eq. [3] reduces to Eq. [5] as expected and the spin signal 
decays exponentially with the distance. Reducing the 
value of the relaxation parameter (of the contact resis- 
tance for given Rsq and W) the length dependence of the 
spin signal starts to deviate from the exponential form 
|32 |. In the case of finite Rc and very long A {R/X 0) 
and short devices {L <^ A), the non-local resistance is 
±{P^RsqR'^)/{WL). Now, the spin signal is indepen- 
dent of A, it depends on the injector detector spacing as 
1/L. Due to strong contact induced spin relaxation, spin 
transport is relevant only on the device length scale, the 
electrodes effectively cut the graphene flake with respect 
to the spin transport. 



a) L = 5 urn 




L = 5 /im and two values for A, 2 and 10 ^m, respec- 
tively. The signal is normalized to the value correspond- 
ing to infinite contact resistance. For given graphene 
properties, a variable R parameter represents a variable 
contact resistance Rc- The full circles in Fig. [5^) corre- 
spond to i? = 10^^ m (i.e., Rc — 1 kil for typical values 
Rsq = 1 kfi and W — 1 fim). Due to the conductivity 
mismatch, for A = 2 /xm one measures only 25 % of the 
totaly available signal in the ideal case of high impedance 
contacts. For A = 10 ^m, only 3.7 % of the ideal signal 
is measured. The strong contact induced spin relaxation 
regime is identified as the region where the signal is pro- 
portional to R^. In Fig. [Sla) we show the variation of 
the spin signal with respect to the injector-detector sep- 
aration for A = 2/im and three values of the relaxation 
parameter. For R = lO^^m and L < X the spin signal 
depends on the injector-detector separation as 1/L. 



B. Graphene spin valve measurements 

In this section we discuss the general characteristics 
of the spin valve measurements. In Fig. [6] we show spin 
valve measurements for a type I and a type II device for 
three charge density regimes: zero charge density (Dirac 
neutrality point) and metallic hole and electron densi- 
ties of 2-3x 10"'^^cm~2. Due to the different contact sep- 
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FIG. 5: The variation of the spin valve signal (normalized to 
the maximum spin signal ideally available) as a function of: a) 
the relaxation parameter R for two values of A = 2, 10 /im; b) 
injector-detector separation L for R = 10"''^, 10~^, 10~^ m. 

In Fig. [5^) we plotted the variation of the non-local 
spin valve signal with respect to the relaxation param- 
eter R for a device with an injector-detector separation 



FIG. 6: Spin valve measurements for a type I device (left col- 
umn) and a type II device (right column), dimensions given 
in the insets in micrometers. For the Dirac point (middle 
graphs) the spin signals are smaller than for metallic regimes 
(densities of 2-3x lO^'^cm"^). The y-scales have the same span 
for each device. The resistance levels Rl, R2, R3 and R4 ob- 
served for the device II. B are due to the magnetic configura- 
tion of all four electrodes (see Fig. [7] and the text). 
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arations, carrier mobilities and metallic charge densities, 
a direct comparison between these measurements is not 
trivial. In general, we observe large signals (few ohms 
up to 50 ohms) for the type I devices, whereas for the 
type II the spin signals were a few tenths of an ohm up 
to a few ohms. We identify two reasons for this differ- 
ence. First, for the type I devices the contact resistances 
were in the order of 10 to 200 kil, whereas for the type II 
the contact resistances were always lower, of the order of 
1 to 10 kr2. As already discussed the low impedance con- 
tacts induce significant spin relaxation, which severely 
reduces the magnitude of the spin signal. Second, the 
measurements are consistent with relatively high spin in- 
jection/detection efficiencies (P=:10-20%) in the case of 
high impedance contacts, whereas for low impedance con- 
tacts the efficiencies were of the order of 3-5%. 




FIG. 7: Schematic representation of the electrochemical po- 
tentials for the spin up and spin down channels as a function 
distance for the device II. B in Fig. |6l The magnetic config- 
urations of all four ferromagnetic electrodes (which give the 
multiple resistance levels Rl, R2, R3 and R4) and the elec- 
trode positions are indicated by the arrows. 



For the type II device in Fig. [5] we observe four re- 
sistance steps which we associate with the magnetiza- 
tion switching of all four electrodes, graphically shown 
in Fig. [T] Let us consider that for high positive mag- 
netic field the magnetization of all electrodes points in 
the "up" direction and that F2 injects spin up carriers in 
graphene. At Fl, where the current flows from graphene 
to Fl, spin down injection takes place. This partially can- 
cels the spin up accumulation created by F2. Both F3 
and F4 electrodes probe the spin up channel and a posi- 
tive non-local voltage (resistance) is measured since F3 is 
closer to the injectors than F4, Fig. [7^). Ramping down 
the field towards -50 mT the electrodes switch their mag- 
netization in the reverse order of their width: Fl, F4, F3 
and F2. When Fl changes its magnetization orientation, 
both Fl and F2 inject spin up carriers and the total spin 
accumulation increases causing a positive change in Rni , 
Fig. [71d). Next, the detector F4 switches, probing the 
spin down channel, and the non-local voltage difference 
increases again, Fig.[71:). At around -25 mT, F3 switches 
its magnetization. Now, both injectors inject spin up 
carriers whereas the detectors probe the spin down ones 
which causes Rni to become negative, Fig.[7ji). The last 
switch is that of F2. This configuration (not shown) is 
equivalent to the starting one with all magnetizations 
pointing down now and the starting level of non-local 



resistance is recovered. 

The four changes in the non-local resistance point out 
that the spin relaxation length in graphene is at least of 
the order of the F1-F3 or F4-F2 distances (~ 2.9 fim for 
this sample) since otherwise the resistance step due to 
Fl or F4 switching their magnetization would not be ob- 
served. Additionally, the spin accumulation can extend 
underneath all contacts, which means that the aluminum 
oxide barriers are not too opaque nor too transparent. In 
the next sections we extract quantitative information re- 
garding the spin relaxation length A. 



C. Length dependence of the spin signal 

Eq. [3] can be used to extract A from the length depen- 
dence of the spin signal. Ideally, all one has to do is to 
make a series of devices which should strictly meet the 
following conditions. First, the Fl and F4 electrodes do 
not contribute to the measurement, they are either placed 
at distances much larger than A with respect to F2 and 
F3 or they are replaced with non-magnetic ones. Second, 
F2 and F3 should have identical spin inject ion/ detect ion 
efficiencies and resistances from device to device. In 
addition, all the graphene flakes should have identical 
transport properties (same carrier mobility) and the same 
width (in case that enhanced spin relaxation takes place 
at the device edges). The samples should be processed 
and measured in identical conditions, since it is not ac- 
tually known how the substrate could influence the spin 
relaxation. If all these conditions are met, then the only 
variable left is the spacing between the injector and the 
detector electrodes. This has to be varied systematically, 
so that the application of Eq. [3] would allow one to extract 
A and P. Experimentally some of these conditions can be 
fairly met (especially if several devices can be fitted on 
the same flake) , except controlling the polarization of the 
electrodes and the reproducibility of contact resistance 
which are probably the most important parameters. In 
spite of these problems we have been able to show some 
consistent behavior for sarnples where the contact resis- 
tances were more uniform Q, see also Fig. [TUk). 

In this paragraph we show how the multiple resistance 
steps, which are due to ferromagnetic electrodes placed 
at various distances, can be used for a quantitative esti- 
mation of A. For the all up magnetization conflguration 
which gives the resistance level Ri (Fig. [6] and Fig. [7^) 
the spin accumulation present at the position of the volt- 
age probe F3 is just the sum of two spin accumulations of 
opposite sign produced by F2 and Fl which have decayed 
exponentially with the distance (F2-F3 and F1-F3). We 
neglect here the influence of the electrode F2 on the spin 
accumulation produced by Fl and detected by F3. The 
same arguments apply for the detector F4 but the expo- 
nential decays are for larger distances. Additionally, we 
note that some measurements (including this one) show 
a gate voltage dependent background resistance (which 
we denote as Rq). Note that, theoretically, the spin sig- 
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FIG. 8: Length dependence of the spin signal as extracted 
from the magnitude of the multiple resistance levels (type II 
devices) for three conduction regimes: n~ and hole/electron 
densities of 2-3x lO^'^cm"'^. For these samples we calculated 
an average R parameter of l-2xl0~^ m. The spin relaxation 
lengths are higher in the metallic regime (subscript e/h for 
electron/hole conduction) than at the Dirac neutrality point 
(subscript D). 



nals i?2 and R4 should have identical values but oppo- 
site signs, see Fig. [7]b) and d). Therefore, we can write 
Ri = [+R23 - ^13] " [+R2i - Rii] + Ro, where R23 
is the resistance level produced by the F2-F3 injector- 
detector pair which are placed at the distance L32 and 
so on. Applying a similar reasoning for R2,R3 and i?4 
we obtain a system of four equations with five unknowns 
R23, Ri3,R24, Ri4 and Rq. In this sample the F1-F3 and 
F2-F4 distances were virtually equal (~ 2.9 /im) so by 
taking i?i3 — R24 we can solve the system of equations 
and extract the length dependence of the spin signal from 
a single measurement. The same analysis has been per- 
formed for another sample. The results are plotted in 
Fig. [H] together with fits according to Eq. O The fits 
reveal a spin relaxation length of the order of 1 /im, 
somewhat higher in the metallic regime than for the zero 
charge density regime. 



D. Modeling the contact resistance effects in the 
spin precession measurements 

In this section we discuss the effect of the contact 
impedance on the spin precession measurements. Let 
us consider a spin valve device in which the ferromag- 
netic injector/detector electrodes have their respective 
magnetizations oriented parallel and a magnetic field , 
smaller than what is required to pull the magnetization 
of ferromagnetic electrodes out of plane, is applied per- 
pendicular to the sample plane. In this situation, the in- 
jected spins undergo precession (Hanle effect) while dif- 
fusing from the injector towards the detector, see the 
schematics of the inset in Fig [3^) . The detector signal 
oscillates with B^- The signal is maximum at Bz — 
(no precession), vanishes at fields which correspond to 
an average spin precession angle of 90 degrees and it is 
minimum when the spins have precessed an average angle 
of 180°. In Fig. [S] we present several calculated preces- 
sion signals for various values of the relaxation parameter 



a) 




D = 2e-2 mVs 
T = 200 ps 
X = 2 ^m 



-R= 1e+9 m 
-R=1e-6 
-R= 1e-9 



b) ■ 


^ 1 1 

R = 1 e-9 m . 




A = 2 urn 




D= 


0.5e-2 mVs, t= 800 ps 


D= 


2.0e-2 T= 200 


D= 


8.0e-2 T= 50 



c) 




R = 1e-9 m 
T = 200 ps 



0.5e-2 m /s 

2.0e-2 

8.0e-2 



d) 


— 1 ' 1 

R = 1e-9m - 




D = 2e-2 m^/s 




T= 50 ps 




T= 200 




T=800 



0.4 



B(T) 



0.2 
B(T) 



FIG. 9: Calculations of the spin precession curves for various 
values of R, D and A. Only the parallel magnetization of the 
injector/detector electrodes case is shown for positive values 
of the perpendicular magnetic field. 



R, the diffusion constant D and the spin relaxation time 
T. All relevant parameters are indicated on the figures. 
To allow comparison the signals have been normalized to 
the maximum value. We refer to the Appendix for mod- 
eling details. Next, we analyze the shape of the curves 
with respect to their width, defined by the 90° precession 
angle, and the amount of overshoot, defined as the ratio 
between the minimum and the maximum values of the 
signal (180° and 0° precession). In Fig. [5^) we plotted 
the precession signal for three values of R and for fixed 
D and r (values close to the ones corresponding to our 
samples). With decreasing the relaxation parameter i?, 
i.e. decreasing the contact resistance for given graphene 
properties, the precession curve becomes broader and the 
overshoot increases in magnitude. Of interest is to deter- 
mine how i?, D, r influences the shape of the precession 
curve, and whether one can accurately extract the spin 
relaxation time in the case of very low impedance con- 
tacts, for example Rc — l^i, which for Rsq ~ Ikfl and 
W — 1/im corresponds to R = lO^^m. 

In Fig.[3]D) we show three curves for which R — lO^^m 
and the values for D and r are such that they give a con- 
stant spin relaxation length of 2 /im. This plot indicates 
that the overshoot depends only on A. Extensive anal- 
ysis reveals that for a given injector-detector distance L 
and for X < L, the overshoot depends only on A and R, 
it does not depend on the individual values of D and r 
as long as they produce the same A. Therefore, for a 
given sample (for which L and R are known) the spin re- 
laxation length can be estimated just by quantifying the 
amount of the overshoot. In Figs.[5J:) and d) we plotted 
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spin precession signals by varying either D or r while 
keeping the other parameters constant. For a given r, 
the precession curve is wider for a higher D. This can be 
simply understood by the fact that for a higher D (i.e. a 
higher carrier mobility), a stronger magnetic field is re- 
quired to induce 90° precession since it takes less time for 
the spins to diffuse the fixed injector/detector distance. 
Fig. [SJi) shows wider precession curves for smaller spin 
relaxation times. For high relaxation times, the shape of 
the curves is essentially the same. This can be under- 
stood by the fact that most of the injected spins do not 
relax before reaching the detector and the angle of pre- 
cession is mainly determined by the time it takes them 
to arrive at the detector, i.e. the diffusion constant and 
the device length. For low spin relaxation times, the 
spins following a long diffusive path have a high chance 
to relax before reaching the detector and, therefore, do 
not contribute to the measured signal. In this case, the 
spin signal is mainly determined by the spins following a 
short diffusive path, and which arrive faster at the detec- 
tor (effectively having a higher diffusion constant). The 
overall behavior is that, for given device characteristics 
and dimensions, the precession measurements allow accu- 
rate extraction of the spin transport properties provided 
that the contacts do not induce strong spin relaxation, 
i.e. the relaxation parameter R is comparable or higher 
than the spin relaxation length. 



E. Spin precession experiments 

In Fig. [TOb ) we show the length dependence of the 
spin valve signal for three type I devices made on the 
same graphene flake measured at 4.2 K in the metal- 
lic hole conduction regime n{h) ~ 2 x 10^^ cm~^. For 
this sample the contact resistances were more uniform 
(1-2 kil). The measurements indicate a spin relaxation 
length of 1.6 fxin. In Fig. llOb ) we present a spin pre- 
cession measurement for a 5 /um long device from the 
same batch of samples. The fit of the low B field part of 
the precession curve indicates a spin relaxation length of 
1.4 /im, in very close agreement with the length depen- 
dence measurement. This indicates that the longitudinal 
and transversal relaxation times (Tl and T2) are simi- 
lar for these devices [1]. Therefore, we argue that the 
precession measurements can be used to determine the 
spin transport properties (r,A) of graphene. The advan- 
tage over the difficult to realize length dependence of the 
spin signal (due to the irreproducibility of the contact 
resistances or the spin injection/detection efficiencies) is 
obvious: the fabrication and proper characterization of a 
single device. 

With increasing the perpendicular magnetic field above 
0.5 T the magnetization of the Co electrodes is slowly 
pulled out of plane. For fields above 1.25 T the rela- 
tive orientation of the magnetization of the Co electrodes 
is parallel (both magnetizations are completely out of 
plane). Now, the injected/detected spins are perpendic- 
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FIG. 10: Measurements at 4.2 K in the metallic hole regime, 
type I devices, a) length dependence of the spin valve signal 
and b) spin precession. 



ular to the two dimensional plane of the graphene layer 
and parallel to the external field (no precession takes 
place) . We observe that the spin signal is somewhat lower 
than the zero field signal, when the orientation of the in- 
jected/detected spins is parallel with the sample plane. 
The difference in the magnitude of the signal corresponds 
to a relaxation time t± smaller by ~ 20 % than r|| . 
Similar room temperature experiments not shown here 
support these low temperature spin anisotropy measure- 
ments. This indicates that the effective magnetic fields 
responsible for the spin relaxation lie mostly in the two 
dimensional plane of graphene. For a full discussion we 
refer to |8|. However, we note that in the case of Rashba 
or Dresselhaus type spin orbit interaction the effective 
spin orbit fields lie exclusively in the graphene plane and, 
theoretical calculations indicate that t± = (1/2)t|| [sl]. 

In Fig. [TT] we show spin precession measurements for 
a type I (contact resistances of 20-40 kfl) and a type II 
sample (contact resistances of 1-2 kfi) for similar charge 
carrier type and density. Even though the mobilities and 
widths of the two devices were different, the spin relax- 
ation lengths are rather similar. This indicates that the 
presence of the aluminium oxide layer on top of graphene 
(type I sample) does not influence the spin transport 
properties signiflcantly. 

Overall the values of the diffusion constant agree well 
with the ones extracted from the conductivity measure- 
ments (see [3, 8] for details), whereas the spin relaxation 
times are of the order of 100-200 ps. Only a few devices 
show relaxation times smaller than 100 ps. In all the mea- 
surements we observe that the spin transport properties 
depend on the gate voltage, i.e. on the charge density. 
Generally, for the same device we observe that both the 
spin diffusion constant D and the spin relaxation time 
T at the Dirac point are approximately a factor of two 
smaller than at high carrier densities, see Fig. 1111 Since 
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FIG. 11: Room temperature spin precession measurements 
for a type I device (left column) with _Rc=20-40 kQ, fi — 
2.4 X 10^ cm^/Vs and for a type 11 device (right column) 
with Rc=l-2 kf2, /i = 5 x 10^ cm'^/Vs in different conduction 
regimes. The carrier concentrations are expressed in units of 
10^^ cm~'^. The flake dimensions are indicated on the figure. 
The solid red lines through the data points are fits with the 
parameters indicated on the figures. D, r and A are expressed 
in units of 10~^m^/s, ps and fim respectively. 



D = VFTe/2 (with vp the Fermi velocity and the mo- 
mentum relaxation time), it seems that r oc Tg. This 
behavior is expected for the Elliott- Yafet type spin scat- 
tering mechanism. 
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FIG. 12: Room temperature spin precession measurements in 
the metallic hole conduction regime n{h) ~ 2 x 10^^ cm"'^ for 
four type III devices of different widths. 



2.5x10^ cm2/Vs, the spread in the spin relaxation times 
is rather large. Generally, the spin valve signals at the 
Dirac point were too small to produce a useful preces- 
sion curve. For some precession curves there are a few 
sets of parameters that describe reasonably the experi- 
ment. In that case we fixed the spin diffusion constant 
to the charge diffusion constant, leaving the spin relax- 
ation time as the only relevant fitting parameter. In the 
order of increasing the graphene flake width (0.10, 0.24, 
0.39 and 2 /im) the spin relaxation times read 40, 200, 
154 and 60 ps respectively. No clear scaling behavior is 
found. Though the spin relaxation time for the narrowest 
device (100 nm) shows the lowest value of 40 ps, a defini- 
tive conclusion with respect to a presumably strong spin 
relaxation at the device edges cannot be drawn since spin 
relaxation times of 50-60 ps were measured for relatively 
wide flakes as well. 



F. Spin precession in graphene ribbons 



IV. CONCLUSIONS 



In this section we discuss room temperature spin pre- 
cession measurements for structures down to 100 nm in 
width (type III devices). These measurements were per- 
formed in order to establish whether the short spin re- 
laxation times may be due to enhanced spin flip pro- 
cesses taking place at the edges of the graphene flake. 
The fabrication procedure involved an oxygen plasma 
etching step and it was described in the Experimental 
section. The values of the contact resistances for this 
type of etched structures show the poor reproducibil- 
ity compared with the type I and II devices. We at- 
tribute this to the unknown residues produced during the 
etching process which hinders making proper contacts to 
graphene. In Fig.[T2]we show a set of four precession mea- 
surements performed on different devices in the metallic 
hole conduction regime. In spite of similar charge trans- 
port properties, with carrier mobilities of the order to 2- 



We have successfully realized all electrical injection and 
detection of spin accumulation in graphene at room tem- 
perature and 4.2 K. The conductivity mismatch prob- 
lem has been partially overcome by the introduction of 
a thin aluminium oxide layer between the metallic Co 
electrodes and the semiconducting graphene. By proper 
modeling, we have shown that the spin relaxation via the 
contacts does not dominate the spin transport properties 
of our devices. From the length dependence of the spin 
valve signal and spin precession measurements, we have 
shown that the longitudinal and transversal spin relax- 
ation times are similar. For the 2-5x10"^ cm^/Vs carrier 
mobilities of our samples, we found spin relaxation times 
of 50-200 ps. Good agreement has been found between 
the spin diffusion constant and charge diffusion constant. 
The linear dependence of the spin relaxation time on the 
momentum scattering time indicates that the spin relax- 
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ation mechanism is Elliot- Yafet type. For carrier mo- 
bilities a factor of 100 higher than our samples, a spin 
relaxation length of the order of 100 micrometers is ex- 
pected. The anisotropy in spin relaxation times reveal 
that the effective magnetic fields due to the spin orbit 
interaction lie mostly in the two dimensional graphene 
plane. In the range of carrier mobilities of our samples, 
the scattering at the devices edges, for structures down 
to 100 nm in width, and the presence of the aluminum 
oxide covering the graphene flakes in some experiments 
do not influence the relaxation times significantly. 
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Appendix. Modeling of spin relaxation induced by 
the contacts 



We consider the geometry shown in Fig.[T3]with Fl and 
F4 being non-magnetic. The injector F2 and the detector 
F3 are placed at the distance L. In the frame of the two 
channel current model the current is carried by two inde- 
pendent channels, a spin up and a spin down one. In the 
absence of external magnetic fields the spin accumulation 
is constant in the y and z directions, only the x com- 
ponent varies with the distance. To simplify the math- 
ematical description we consider symmetric splitting of 
the electrochemical potentials of the two spin species and 
discuss only the spin up channel. Neglecting the linear 



F2 F3 



SiO, 



n-doped Si 



x=0 



x=L 




FIG. 13: Modeling of the spin relaxation via the contacts. Fl 
and F4 are considered as non-magnetic. 

term due to charge current flow for x < 0, the solutions 
for the spin up electrochemical potential in the x direc- 



tion are of the form 



li-\{x) = aexp(-|- — ), for a; < 
A 

X X 

Hl{x) — 6exp(— — ) -f cexp(-|- — ), for < a; < L 
A A 

X 

Ulix) = dexp(— — ), for x > L, 
A 

Considering parallel orientation of the magnetization of 
the injector/detector electrodes and the same polariza- 
tion P for the injector/detector electrodes, we write the 
non-local resistance as 



Rnl {L) 



P _ P deyi^{~L/X) 



(6) 



where we took into account the symmetric splitting 
fJ.]{L) = —fj,i{L) when applying Eq. 2 from [3^. The 
constants a, b, c, d are determined from the boundary con- 
ditions. The continuity of the spin up electrochemical 
potential at x ^ 0, L gives the following two equations 



(7a) 



L. 



dexpi-j). (7b) 



The spin relaxation via the contacts is taken into ac- 
count when writing the spin current conservation equa- 
tions as explained in the following. The ferromagnet F2 
injects a spin polarized current in graphene via tunnel 
barriers of resistance Rc with an efficiency P. The to- 
tal injected spin current is ls{0) = Plac and the total 
spin up current available is half, /|(0) — Placf^ = lx{0)- 
The spin up current lx{0) diffuses to the left and to the 
right in graphene or flows back into the contact. The 
spin up current flowing in graphene is written in the 
form Iq = ±{a/2)W{l/e)dii^/dx. Here e is the elec- 
tron charge, a and W are the conductivity and the width 
of the graphene fiake. The (unwanted) spin up current 
through the contact is written as ij-, = (l/i?J,)^|/e « 
(2ei?c)^^MT- Here we have considered R[;{— Rq) ~ 
2Rc, with i?p the effective resistance the up spins en- 
counter when going back into contact. The spin current 
conservation equations at a; = 0, L read 



plac_ 



a 



aW 
2^ 



aW 



aW 



2eA 2eA 2ei?c 



(8a) 



= 



,crVF , aW , L, 

^2a^"P(-A) + ^2a^"P(+A) 



"^2a^"P(-A)+2ei?c 



d , L , , , , 

exp(--) (8b) 



Here we have used the fact that the spin current flowing 
away from the contact (the exponential decaying elec- 
trochemical potential from the contact point of view) is 
positive. The boundary conditions give a system of four 
equations with four unknowns from which the constants 
a,b,c,d are extracted. Using a = l/Rgq, we find the 
expression for Rni (L) , given by Eq. [31 
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In the presence of an external magnetic field applied 
in the z direction the general solutions of the spin up 
electrochemical potential are of the form [13] 



0,L 



J^J- I I cxp(±7^ix) 







(9) 



0,L 



1 

—i ] exp{±K2 x) 




C'^^ I i I exp(±if+a;) 



where the — (+) signs correspond to the positive (neg- 
ative) X direction and Ki and are given by Ki = 
(%/I>r)-i = andK^ = A~Vl ± «wr. Since the elec- 
trochemical potential cannot increase when x approaches 
— (-t-)oo, the solutions for x < {x > L) contain only the 
exponentially decaying terms, i.e. the terms containing 
+Ki,2 for a; < and the terms —Ki,2 for x > L. For 
< a; < L the solutions contain six terms so that the 
back flow of spins due to the presence of the detector 
is taken into account. There are in total twelve coef- 
ficients {Al,Bl,Cl, A°_,B°_,C°_,A^,B'^,C^,A^,B^,C^) 



that describe the solutions. These are found from the 
boundary conditions: the continuity of electrochemical 
potentials and the spin current conservation, laws which 
we apply for each x,y and z component. The continuity 
of the electrochemical potentials give in total six equa- 
tions, three for injection and three for detection. The spin 
current conservation laws give another six equations and 
they are written in a similar manner as it was described 
above. Except the x component lx{0) all the other spin 
current components ly{0), lz{0), Ix{L), Iy{L), Iz{L) are 
null. The spin up electrochemical potential present at 
the detector (x component) is given by 

/i|(L) B^ exp i-K^L) + exp {-K+L). (10) 



The boundary conditions give a system of twelve equa- 
tions with twelve unknowns. For given r and D we solve 
the system of equations numerically to find the coeffi- 
cients i?:^ and Ch for each value of the magnetic field. 
P is just a multiplication factor. We calculate a series 
of precession curves for a mesh given by different values 
of r and D. From this set of curves we choose the one 
that shows the smallest deviation from the experiment by 
the least mean squares method. The parameters r, D, P 
corresponding to this curve represent the best fit. 
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